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ABSTRACT 

We report the first results of a spectroscopic search for Lya envelopes around three z ~ 4.5 radio-quiet quasars. Our observational 
strategy uses the FORS2 spectrograph attached to the UT1 of the Very Large Telescope (VLT) in the multi-slit mode. This allows us 
to observe simultaneously the quasars and several PSF stars. The spectra of the latter are used to remove the point-like quasar from 
the data, and to unveil the faint underlying Lya envelopes associated with the quasars with unprecedented depth. We clearly detect 
an envelope around two of the three quasars. These envelopes measure respectively 10" and 13" in extent (i.e. 67 kpc and 87 kpc). 
This is 5 to 10 times larger than predicted by the models of Haiman & Rees (2001) and up to 100 times fainter. Our observations 
better agree with models involing a clumpy envelope as in Alam & Miralda-Escude (2002) or Chelouche et al. (2008). We find that 
the brighter quasars also have the brighter envelopes but that the extend of the envelopes does not depend on the quasar luminosity. 
Although our results are based on only two objects with a detected Lya envelope, the quality of the spatial deblending of the spectra 
lends considerable hope to estimate the luminosity function and surface brightness profiles of high redshift Lya envelopes down to 
F rsj 2 — 3 x 10 -21 ergs -1 cm -2 A -1 . We find that the best strategy to carry out such a project is to obtain both narrow-band images 
and deep slit- spectra. 
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1. Introduction 

Large-scale Lya emission (^10-10 kpc) is common a mong 
high redshift (z > 2) radio galaxies ( Vill ar-Martinll2007"|) . This 
extended gas shows two components: one is kinematically per- 
turbed by the jets and is part of a jet-induced outflow; the other 
has a lower velocity dispersion (a few hundreds of kms -1 in- 
stead of about 1000 km s -1 ) as well as a fainter surface bright- 
ness, but may extend beyond 100 kpc. It has been shown, at least 
in the case of M RC 2104-242, that t he gas is infalling towards 
the galaxy center dVillar-M artm 2007). This author makes the in- 
teresting suggestion that "the radio activity is fed by the infalling 
gas, so that it is only detected when the infall is happening and 
efficiently feeding the active nucleus". 

Similar Lya envelopes (sometimes called "blobs") are 
also found around rad i o-quiet quasar s (ISteidel et al.l 12000: 
van Breugel et al.1 120061: iBunker et al.1 120031: IWeidiger et al.1 
2005; Christense n et al.l I2006L hereafter CJW). According to 
the latter authors, Lya envelopes around radio-quiet quasars 
(RQQs) are an order of magnitude less luminous than those 
around RLQs, presumably because the emission of RLQ gaseous 
en velopes is enhanc e d by i nteractions with the radio jets. 

Haiman & Rees ( 2001) have envisaged Lya emission as a 
possible constraint on galaxy formation. Gas infalling into the 
gravitational well of a dark matter halo would be heated, and dis- 
sipate this thermal energy partly by collisional excitation of the 



Table 1. Journal of observations, along with the main charac- 
teristics of the quasars. The apparent magnitudes are given in 
the AB system. They are computed by integrating the quasar's 
spectrum through the RSPECIAL ESO filter curve. The absolute 
magnitude assumes Ho = 72 km s _1 Mpc _1 and (Q m ,Q\) = 
(0.3,0.7). 



* Based on observations made with the FORS2 multi-object spectro- 
graph mounted on the Antu VLT telescope at ESO-Paranal Observatory 
(programme 079.B-0132B; PI: P. North) 
** Chandra Fellow 
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Lya transition. But the Lya emission would remain too faint 
to be detectable at high redshift. However, if a central quasar 
turns on and photoionizes the gas, the Lya emission will be 
much enhanced and should be detectable with present-day in- 
struments. For quasars at redshift z = 3 — 8, these authors 
predict surface brightnesses in the range 10 -18 < \i < 10 -16 
ergs -1 cm~ 2 arcsec -2 and ang ular sizes of a few arcseconds. 
Alam & Miralda- Escudd ([2002) have proposed another theoret- 
ical prediction of the Lya surface brightness of a quasar host 
galaxy at z = 3. They find /i = 10~ 17 - 5 and fi = 10~ 19 - 5 
ergs -1 cm -2 arcsec -2 at angular separations of 0.5" and 3" re- 
spectively from the quasar. This is mu ch more pessimistic than 
the estimate by Hai man & Reesl ([2001), because the authors as- 
sume a smaller clumping factor of the cold gas and a Lya emit- 
ting region much smal ler than the virial radius. More recently, 
Chelouc he et al.l (120071) suggested that quasar nebulae are an- 
other manifestation of metal absorption systems associated with 
L* galaxies which, by virtue of a nearby quasar, become more 
efficient Lya emitter that also scatter Lya photons from the 
broad line region of the quasar (BLR). The extent and luminos- 
ity of the halo in their model serves as a means to study the 
nearby environments of quasars and weigh the gaseous content 
of their halos. To summarize, the study of the emission prop- 
erties of quasar envelopes provides a promising new means for 
testing galaxy formation models and various scenarios for the 
enrichment of the inter-galactic medium. 

To explore in more detail the spatial extent, the luminos- 
ity and kinematics of the large hydrogen envelopes of remote 
quasars, we have selected a sample of quasars at redshift z ~ 
4.5, spanning a broad magnitude range. Before carrying a sys- 
tematic study of the Lya envelopes of these quasars, we have 
obtained deep VLT optical spectra for 3 of them, spanning 3 
magnitudes. The present paper describes our observational strat- 
egy and our main results, with the detection of a Lya "nebula" 
for two of the three quasars. 
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Fig. 1. Spectrum of the sky emission, extracted from our FORS2 
data and grism G1200R in combination with a 2 // -slit. The full 
wavelength range accessible to this grism is displayed in the up- 
per panel. A zoom is displayed in the bottom panel, of the region 
where falls the Lya emission redshifted at z~4.5. The positions 
of the Lya line of the three quasars are also indicated. 
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2. Deep VLT optical spectroscopy 

Motivated by the hypothesis that the Lya emission around 
quasars is enhanced by the quasar radiation field, our goal is 
to characterize the properties of the Lya envelopes as function 
of the quasar luminosity. The construction of a sample is a com- 
promise between this scientific goal and the technical constraints 
imposed by the intrinsic faintness of the Lya envelopes. The ob- 
servations presented here are intended to demonstrate the feasi- 
bility of the project and to give a first characterization of the Lya 
envelopes for three quasars with different luminosities. 

2. 1. Sample and observations 

We select bright quasa rs from the 3rd editio n of the SDSS 
catalogue of quasars (iSchneider et al.l I2005I) and from the 
12th edition of the Verons' catalogue (iVeron-Cetty & Veronl 
120061) . All our targets are at high redshift, i.e., when 
galaxy formation was still in place. SDSS J0939+0039 and 
Q 213 9—4324 have not been ex plicitly related to any radio 
source (iVeron-Cetty & Veronl 120061) but given their brightness, 
they are likely RQQs BR 10 33- 0327 is an RQQ acco rding to 
IVeron-Cettv & Veronl (l2006h and lBechtold et al.l d2003b . 

A key observational issue is to avoid both the telluric ab- 
sorption lines and the numerous atmospheric emission lines in 
the red part of the optical domain. Residuals of such lines in the 
sky- subtracted spectra may mimic a Lya envelope. This latter 
consideration restricts the accessible redshift range to 4.44 < 
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Fig. 2. The combined and sky- subtracted FORS2 spectra for the 
three quasars, prior to any deconvolution (see text). The images 
have been binned by 5 pixels in the spectral direction (i.e., the 
new pixel size is 3.8 A) and by 2 pixels in the spatial direction, 
resulting in a pixel size of 0.5". The extended Lya emission is 
indicated with arrows in two of the objects. 



z < 4.61, as shown in Fig.Q] with an additional small window 
at 4.26 < z < 4.30. Finally, there must be adequate nearby stars 
with about the same brightness as the quasar in order to carry out 
a careful subtraction of the quasar's l ight u sing image deconvo- 
lution, as described in lCourbin et al.l (l200Qh . 

The observations were carried out in ESO Period 79 (April- 
September 2007) in service mode with the FORS2 multi-object 
spectrograph attached to VLT-UT1 . The ESO grism G1200R+93 
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Fig. 3. The three combined and sky- subtracted FORS2 spectra along with their spatial deconvolution. The spatial resolution after 
deconvolution is 0.25". In each panel, the vertical arrow indicates the position of the Lya emission line of the quasar. A Lya 
envelope is detected in two quasars out of three. The envelope of BR 1033—0327 is redshifted compared with the quasar. The height 
of the spectra is 16" in all panels. 



has a resolving power R= 1070 with a 2" -slit, ensuring to catch 
most of the flux of the Lya envelope. This grism is used in com- 
bination with the GG435+81 order separating filter, leading to 
the wavelength coverage 6000 A < A < 7200 A. The maximum 
efficiency of this combination coincides well with the expected 
wavelength of the redshifted Lya line, i.e. about 6686 A. 

The multi-slit MXU mode is used, with slits that are long 
enough (typical length: ~ 20") to reliably model and subtract 
the sky emission. Although only one scientific target is observ- 
able in each field, we use the MXU capability in order to ob- 
serve several stars through identical slits. In this way, a spectrum 
PSF is measured simultaneously with the quasar. This is cru- 
cial for the spatial de convolution of the data to work efficiently 
dCourbin et al. 2000) and to separate well the quasar spectrum 
from that of the putative envelope. 



Although the requested maximum seeing is 0.9" for our ob- 
serving programme, this condition is fulfilled in practice for 
three Observing Blocks (OBs) out of the nine executed. Each 
planned OB consists of one short acquisition image, followed by 
two long exposures for spectroscopy. The planned total exposure 
time for each object is 2.7 hours, with seeing better than 0.9". 
The journal of the actual observations is presented in Table 1. 

2.2. Reduction and spatial deconvolution 

The data reduction is straightforward. It is carried out using 
the standard I RAF procedures. The individual spectra listed in 
Table 1 are flatfielded using dome flats, and then wavelength- 
calibrated in two dimensions in order to correct the sky emission 
lines for slit curvature. The scale of the reduced data is 0.76 A 
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Fig. 4. Extracted ID spectra for the three quasars. In each case, the top panel shows the quasar alone and the bottom panel shows 
the Lya envelope alone, after spatial deconvolution of the spectra. The vertical dotted line indicates the position of the un-absorbed 
Lya emission line, at the redshift of the quasar. In all cases, the spectrum of the Lya envelope has been smoothed using a boxcar of 
8 A. In the case of BR 1033—0327 we also display the un-smoothed spectrum (middle panel), as the signal-to-noise is higher than 
in the other two objects. 

Table 2. Main properties of the Lya envelopes, as measured on the spectra of Fig. HI All parameters are given in the observer's 
frame. The 1 — a detection limit is the standard deviation of the background noise (after smoothing with a car box of 8 A) integrated 
along the whole slit. Note that the surface brightness is integrated in wavelength but given per arcsec 2 , while the 1 — a limit is 
spatially integrated but given per A. 
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[erg s" 1 cm -2 "~ 2 ] 
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50 ±10 


13, 86 


7.7(±0.8) x 10~ iy 
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Table 3. Lya luminosity of the quasar in the BLR, compared 
with the lyminosity of the extended envelope. The flux in the 
quasar BLR is measured in the wavelength interval 1200-1230A 
(rest- frame). The last two lines give the Lya flux of the enve- 
lope, after correction by slit-clipping (see text). 



Object 


L(BLR) 

[ergs -1 ] 


L(Lya) 
[ergs x ] 


BR 1033-0327 
Q 2139-4324 
SDSS J0939+0039 


7.2(±0.4) x 10 45 
4.6(±0.2) x 10 44 
4.1(±0.2) x 10 44 


4.0(±0.4) x 10 4 ^ 
3.0(±0.6) x 10 40 


BR 1033-0327 
Q 2139-4324 




2.0(±0.2) x 10 4a 
1.2(±0.3) x 10 41 



per pixel in the spectral direction and 0.25" in the spatial direc- 
tion. 

The sky emission is then subtracted from the individual 
frames by fitting a second order polynomial along the spatial 
direction. This fit considers only the 10 pixels on each side of 
the slit and the sky at the position of the quasars is interpolated 
using this fit, both on the quasar and on the PSF stars. 

The cosmic rays are re moved using the L. A. Cosmic algo- 
rithm ( van Dokkum 2001). All the frames are checked visually 
after this process in order to make sure that no signal is mistak- 
enly removed from the data. Particular care is paid to the good 
seeing frames, where the central parts of the quasar and of the 
PSF stars must not be misidentified with cosmic rays. 

Even though FORS2 has an atmospheric refraction correc- 
tor, the shape of the spectra along the spectral direction shows 
slight distortions, i.e., the position of the spectrum changes as a 
function of wavelength. These distortions are corrected for and 
the spectra are eventually weighted so that their flux is the same 
at a reference wavelength, before they are combined into a deep 
2D spectrum. We show in Fig. [2] the combined spectra of the 
three quasars, after binning both in the spatial and the spectral 
directions. An extended Lya envelope is already visible up to 
A" away from the quasars in two objects: BR 1033—0327 and 
Q 2139-4324. 

While a Lya envelope is visible on the combined data al- 
ready, measuring its actual flux requires accurate spatial de- 
blending. We carry out this delicate task by spatial ly deconvolv- 
ing th e spectra following the method described in lCourbin et al] 
(2000). This method is an adaptation to spectroscopy of th e 
"MCS" image deconvolution algorithm (Ma gain et al 1 1 19981) . 
It has been used successfully on many occasions. Among 
others, the algorithm has been used in order to unveil the 
spectrum of the lensi ng galaxy in multiply imaged quasars 
(lEigenbrod et al] 12007) and to study low redshift quasar host 
galaxies dLetawe et al.l 20071) . 

The results of the present application to high redshift quasars 
are displayed in Fig. [3j The algorithm uses the spatial informa- 
tion contained in the spectrum of several PSF stars in order to 
sharpen the data in the spatial direction. At the same time, it also 
decomposes the data into a point-source and an extended- source 
channel. The output of the deconvolution procedure consists of 
two individual spectra, one for the quasar and one for its host 
galaxy (or Lya envelope), free of any mutual light contamina- 
tion. It is therefore possible to estimate the luminosity of the 
Lya emission "underneath" the quasar. Subsampling of the data 
is also possible with the MCS algorithm, hence the pixel size in 
Fig. [3] is half that of the original data, i.e., the new pixel size is 
0.125 // . 



3. Results 

The main scientific information in our data is the Lya luminos- 
ity of the envelopes, their angular size, and their mean surface 
brightness. The measurements of the spectra shown in Fig. [Hare 
presented in Tables[2]&[3j To check the flux calibration of the de- 
convolved spectra, we integrate them in the FORS2 RSPECIAL 
filter. This filter is also used to obtain short acquisition images 
prior to the long spectroscopic exposures. These "spectroscopic" 
AB magnitudes are given in Table 1. We check that they are 
compatible with the simple aperture photometry obtained from 
the short images. 

We detect a Lya envelope in 2 out of 3 objects. Our flux 
limit (integrated over the whole object) is very faint, as indi- 
cated in Table [3 up to two orders of magnitude fa i nter than pre- 
vious studies in this field (IChristensen et al.ll2006l : iBremer et al] 
Note, however, that our limit is computed by integrating 
the spectrum over the full extent of the Lya envelope. 

The measurable extent of the envelopes is r ~ 43 kpc for 
BR 1033-0327 and r ~ 33 kpc for Q 2139-4324, measured 
from the quasar's centroid to the faintest visible isophote, i.e., at 
the detection limit as given in Table O 

While the surface brightness of the Lya fuzz is not affected 
by slit losses, the total luminosity is. Assuming that the Lya en- 
velopes are uniform face-on disks with diameters equal to the 
extents quoted in Table [3 we can estimate the amount of flux 
missed by using a slit width of 2" . The observed luminosities of 
the Lya envelopes for BR 1033-0327 and Q 2139-4324 are 
given in Table [3] as well as the luminosities after the correction 
for the slit clipping. 

The velocity of the Lya envelope in Q 2139—4324 is well 
compatible with that of the quasar. However, the extended Lya 
emission in BR 1033-0327 is redshiftedby AA = 26±5 A, i.e., 
A]/ = +1165±225 km s _1 with respect to that of the quasar. 
We have checked our deconvolution using several PSF stars and 
using different smoothing terms (see lCourbin et al.ll2000l) . These 
checks leave the observed shift unchanged. An explanation for 
this line shift may be that the redshift of the quasar is incorrect, 
as it is measured only from the Lya line itself and from the C IV 
and C III] lines. The Lya emission of the quasar is strongly af- 
fected by the Lya forest, and the carbon lines in quasars are 
known to be blueshifted with respect to the actual redshift of 
the quasar. Using infrared observations, where th e quasar [O HI] 
narrow line is accessible at moderate redshifts, iMcIntosh et al] 
(1999) find an average blueshift of 860 km s 1 and 600 km s 1 
of the C IV and C III] broad lines w ith respect to the \Q III] 
line. A similar conclusion is drawn by Richards et al] (120021) us- 
ing SDSS quasar spectra. The velocity shift we observe between 
the quasar and the Lya envelope of BR 1033—0327, although 
large, is still compatible with a biased redshift measurement of 
the quasar. 

3.1. Redshift dependence 

The observed surface brightness of the two Lya envelopes de- 
tected here is fainter by about 1-2 orders of magnitude com- 
pared to the CJW objects. This difference cannot be explained 
solely by the different redshift ranges of the two samples. While 
our sample is at z\ ~ 4.5, most of the objects in CJW are at 
Z2 ~ 3.3 and so the flux ratio by redshift dimming alone would 
bcF(z 2 )/F(z 1 ) = (l+z 2 ) 4 /(l+^i) 4 - 2.7, i.e., much smaller 
than the observed ratio. 

Taken at face value, our sample as a whole finds that higher- z 
objects tend to have smaller surface brightnesses. Nevertheless, 
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Fig. 5. Top: relation between the total luminosity of the Lya en- 
velopes and that of the quasar in the broad Lya line. Our mea- 
surements are compared with CJW. The open symbols represent 
the direct measurements while the filled symbols are corrected 
for slit clipping (see text). The green symbols show our upper 
limit for SDSS J0939+0039, where no Lya envelope is seen. 
In doing this, we assume a size r=33 kpc (i.e., the mean size of 
the two other objects) and we do not correct for slit-clipping. 
Bottom: the radius of the envelope as a function of L(BLR). 
Note the tight trend followed by the envelopes of the bright 
RQQs of CJW, and the different behavior of the envelopes of 
the fainter quasar Q 2139-4324. Note also that BR 1033-0327 
stands above the maximum size of the objects of CJW, likely due 
to our deeper flux limit. 



it is important to note the possibly strong dependence of the en- 
velope properties on the object luminosity (see below). In par- 
ticular, only one object (BR 1033—0327) in our sample overlaps 



with the luminosity range probed by CJW, and its surface bright- 
ness cannot be rejected at the 3 — o level from being drawn 
from the same sample (where o is the standard deviation and 
the t— test performed on the logarithm of the quantities rejects 
the object at the 94% level). The two other objects in our sam- 
ple do have significantly lower surface brightness, at the 99.8% 
confidence according to t— test. This is a lower limit given the 
upper limit on the undetected envelope in SDSS J0939+0039. 
However, the quasar luminosities are also lower by almost an 
order of magnitude relative to the faintest object in CJW. 

One potential concern is whether our observations artificially 
find larger envelopes than CJW just because they are deeper. We 
note, however, that the mean size of the envelopes in our sam- 
ple is ~ 38 kpc compared to 26.4 kpc in CJW's sample, which 
could account for a 50% difference in area but not for an order 
of magnitude effect. 

3.2. Dependence upon the luminosity of the quasar 

While the surface brightness of the Lya envelopes in our small 
sample of 3 objects is much lower than the observations of CJW, 
their total luminosities agree much better. 

We show in Fig.[5]the luminosity-luminosity diagram, com- 
paring the total flux in the Lya envelope as a function of the 
quasar's luminosity in the broad Lya line. After correcting for 
slit-clipping, we find that our points, combined with the ones of 
CJW, follow the linear relation: 

Log[L(Lya)} = 1.8 x Log[L(BLR)] - 39.2 (1) 

This fit considers a mix of objects at redshifts z ~ 3 and 
z ~ 4.5, i.e., we assume no strong redshift evolution of the lu- 
minosity of the envelopes. We do not include in the fit our upper 
limit for SDSS J0939+0039although the data point is shown in 
Fig. [5j The depth of our observations for this object is consis- 
tent with the trend found by using the 5 points of CJW and our 
2 objects with a positive detection of a Lya envelope. Our new 
observations for the 3 objects therefore seem to support the fact 
that brighter quasars also have brighter Lya envelopes, under 
the assumption of negligible redshift evolution. 

The relation between the size of the envelope and the quasar 
luminosity is much less clear. CJW find a trend that brighter 
RQQs also display larger envelopes. In our sample we find no 
such trend, as shown in Fig. [5j The main difference between 
the observations of CJW and ours is depth. With our deeper ob- 
servations, we probe Lya envelopes much further away from 
the quasar than CJW. In addition the small field of view used 
in CJW implies severe clipping of the envelopes, if they extend 
much beyond a few arcsecs. This may also be at the origin of the 
discrepant surface brightnesses between CJW and the present 
work. 

4. Comparison with models and future prospects 

Current models for the emission of quasar envelopes are in their 
infancy and estimates for the observable properties of such ob- 
jects depend on several tunable parameters such as halo mass, 
clumping facto r for the gas, quasar luminosity, metallicity, as 
well as redshi ft (Haiman & Rees 2001; Alam & Miral da-Escudd 
120021: IChdouche et al.l l2007h . This limitation, as well as the size 
of current samples and the parameter range covered, make a 
quantitative comparison with current models unwarranted. In ad- 
dition, our results may be partly contaminated by emission from 
excited inter stellar medium in the host galaxy of the quasar; an 
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effect not accounted for by current models. Nevertheless, several 
qualitative observations may be made. 

The flux level observed from the envelopes of bright quasars 
seem to indicate that such objects inhabit halos which are more 
massive in terms of their gaseous c ontent (and possibly to- 
tal mass) than those of L* galaxies (lAlam & Miralda-Escudd 
2002i lChelouche et al.l2007l) . Further support for th is conclusion 
comes from study of bright quasars at low redshifts fSerber et all 
(2006). The new results presented here for the less luminous 
quasars ind icate that the surface br ightness is lower than the pre- 
dictions of|5 aiman & Reesl (1200 lb for all relevant halo masses 
(see their Fig. 2). Within the framework of their model this sug- 
gests a rapid evolution of those systems already at z > 4.5, i.e. 
the cold gas has been converted to stars (or into a thin disk) al- 
ready at z ~ 4.5. As a consequence, there is not enough gas 
left when the quasar turns on, to reprocess the quasar light into 
detectable Lya emission. 

Our results seem to indicate a much more complex picture 
than is portrayed by current models. To be able to make some 
progress and to quantitatively compare data and models it is nec- 
essary to gain better information concerning: 

- The surface brightness profile of quasar envelopes. This will 
allow to make a better distinction between recombination 
flux and photon scattering models as well as to understand 
the origin of the gas and possible link with absorption-line 
systems associated with gal axies (lAlam & M iralda- Escudd 
120021: IChelouche et alJl2007h . 

- The gas metallicity (via the detection of e.g., NVA1240 
line). This would provide an essential clue to the origin of 
the gas and whether or not it is considerably enriched. It may 
also reveal whether the cooling time- scales and gas densities 
assumed in different models are realistic. 

- Dependence of envelope observables on quasar properties. 
This will allow us to discriminate between the effect of en- 
vironment and that of quasar luminosity on the properties of 
the envelope thereby allowing us to constrain its physics. 

So far, the observational material available to study Lya en- 
velopes of RQQs is extremely scarce, with a mere half-dozen 
objects with different redshifts and observed with different in- 
struments, down to very different depths. 

Our deep slit-spectroscopy observations show that a homo- 
geneous sample of quasars with redshift up to z=4.5 can be built 
in a reasonable amount of telescope time. In addition, our find- 
ing that Lya envelopes can be surprisingly large (r^lO-15") and 
faint, make integral field spectroscopy a poor solution to under- 
take this task, due to poor cosmetics and sensitivity. A more 
viable option is to conduct a two-step programme focusing on 
1- narrow-band imaging of RQQs at different redshift slices in 
order to map surface brightness profiles for the Lya envelopes 
of quasars with a broad range luminosities and 2- to carry out 
deep slit-spectroscopy with well controled slit-clipping, with 
high spectral resolution and with very deep flux limit. The lat- 
ter shall provide us with the additional metallicity and velocity 
information required to constrain models. 
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